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Abstract. When two non-relativistic particles scatter in one dimension, they can 
become entangled. This entanglement process is constrained by the symmetries of 
the scattering system and the boundary conditions on the incoming state. Applying 
these constraints, three different mechanisms of entanglement can be identified: the 
superposition of reflected and transmitted modes, momentum correlations of the 
reflected mode due to inversion of the relative momentum, and momentum correlations 
in the transmitted and reflected modes due to dependence of the scattering amplitude 
on the relative momentum. We consider three standard potentials, the hard core, 
Dirac delta, and double Dirac delta, and show that the relative importance of these 
mechanisms depends on the interaction and on the properties of the incoming wave 
function. We find that even when the momenta distributions of the incoming articles 
are sharply peaked, entanglement due to the momentum correlations generated by 
reflection can be quite large for particles with unequal mass. 
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1. Introduction 

Before two particles scatter, they are in uncorrelated states. If we assume each particle 
can be described by a pure state when the particles are far apart before the 
interaction, then the initial total state of the system is the product of the two one- 
particle states \4> in ) = \4>\) <g) \4>2)- As the particles approach, the state evolves and 
interparticle separability is lost. Separability does not return even if the interaction 
is elastic and after a long time the particles are far removed beyond the interaction 
region. The boundary conditions of scattering are inherently time-asymmetric from 
the perspective of entanglement: practically, entangled particles cannot be easily 
prepared, and even if they could be, they would not approach, interact, and then 
emerge in a separable state. The amount of entanglement generated by scattering is also 
constrained by symmetry. For non-relativistic scattering, free particles are associated to 
projective representations of the Galilean group extended by mass. Galilean invariant 
interactions imply conservation principles, and these shape the mechanisms by which 
the entanglement occurs. 

For simplicity, we consider the generation of entanglement in the non-relativistic 
scattering of two structureless, distinguishable particles in one dimension. Within 
this limited context, we show that the constraints of symmetry and time-asymmetric 
boundary conditions on the the incoming state imply that scattering entanglement 
proceeds by the combination of three mechanisms. The simplest is the entanglement 
between transmission and reflection. This mechanism of entanglement is also the 
coarsest, because as long as the transmitted and reflected modes are orthogonal, 
each particle can effectively be thought of as a two-level system, with each level 
corresponding to a different side of the interaction region. The other two mechanisms 
create entanglement by distorting the wave functions of the modes in a non-separable 
fashion. If particles have different masses, the wave function of the the reflected mode 
is distorted by the transformation that reverses the direction of the relative momentum. 
In other words, because of the conservation of total momentum, momentum correlations 
between the particles are created by reflection since their momentum distributions 
have a finite width. Finally, the scattering amplitudes typically vary with the relative 
momentum and this causes momentum correlations. The wave function of both the 
transmitted and reflected mode in the outgoing state can be inseparably distorted by 
this effect. The contribution of this mechanism to the entanglement of the outgoing 
state diminishes as each particle's wave function becomes sharply peaked about a central 
value, but we will show that the entanglement due to the reflection mechanism does not. 

The study of how entanglement is generated in scattering has interest for a 
variety of reasons. Some leading possibilities for practical implementations of quantum 
information processes, such as ultracold atoms and some solid state devices, are physical 
systems where scattering is central to the dynamics. Also, quantum information 
theory with continuous variables and mixed continuous-discrete variables has many 
open questions, and scattering systems provide a rich structure for exploration of such 
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systems. Finally, scattering is a fundamental method of interaction for systems at all 
quantum scales, and one could hope that entanglement might provide new perspective 
on this basic interaction process. 

Entanglement generation in non-relativistic scattering of structureless, distinguish- 
able particles has been considered previously [UEJEIIUEJEJITIEJE]. Most of these 
previous treatments consider interactions of a single species [21 IH 13 El [9] , and have 
focused on particular special cases. This article provides a unified treatment of all 
these results and provides a framework for further exploration and generalization. The 
general methods applied here were developed in [10] , in which the effect of linear trans- 
formations of observables on the entanglement of a wave function is described. In [10J, 
this technique is used to prove that entanglement with respect to certain sets of ob- 
servables (such as internal-external entanglement, but not interparticle entanglement) 
is conserved in scattering, and this is applied to the mechanism of entanglement due to 
reflection. 

In section 2, this article will show that by combining the results of pU] with the 
constraints of symmetry and the time-asymmetric boundary conditions of scattering, the 
three mechanisms for entanglement can be identified. In section 3, the interplay of these 
mechanisms is explored for three finite-range potentials: hard core, Dirac delta, and 
double Dirac delta. Two important features emerge that contradict some assumptions 
found in the literature. First, we will show that when reflection dominates the scattering 
interaction, even for very narrow momentum distributions, the mechanism of reflection 
distortion can dramatically increase the entanglement. All that matters is the relative 
masses and momentum variances of the particles in the incoming state, not the overall 
scale, and this effect dominates at low energy. Second, we find that although scattering 
resonances cause rapid variations in the scattering amplitudes and therefore the third 
mechanism is relevant, they do not always increase the total entanglement within the 
transmitted and reflected modes. 

2. Entanglement, symmetry, and scattering boundary conditions 

For two structureless, distinguishable, non-relativistic particles in one-dimension, the 
momentum operators of each particle {Pi,/^} form a complete set of commuting 
observables (CSCO). A tensor product structure corresponding to this CSCO is 
Hi <S> Ti.2, where Hi is the single, free-particle Hilbert space. A pure state of the 
system is described by the bi- momentum wave function (fi(pi,P2)- For simplicity, we 
shall restrict considerations so that <fi(pi,P2) £ <S (IR. 2 ) , the Schwartz space on the 
bi-momentum plane (pi,^)- This means that we consider wave functions that are 
smooth, infinitely different iable, and rapidly decreasing at infinity. With this mild and 
physically reasonable restriction, we can employ Reimann integrals and be assured of 
their finiteness. 

The interparticle entanglement of a general pure state \(f>) G TCi <8> H2 can be 
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(1) 



(2) 



For the continuous- variable wave function 0(pi,p 2 )> the purity of the reduced density 
matrix becomes 



The interparticle purity p^ is an entanglement monotone (more purity always mean less 
entanglement) and takes values in the interval (0, 1] as long as the wave function <p(pi, p 2 ) 
is normalized. For continuous variable entanglement, the purity is useful because, unlike 
the entropy of entanglement, one does not need to diagonalize the reduced density matrix 
to calculate pi 2 numerically. Additionally, the simple form ([3]) allows for analytic results 
in certain cases (see below). 

Without loss of generality, all calculations can be performed in the center-of-mass 
(COM) reference frame where the expectation value of the total momentum operator 
P = Pi <g> I 2 + Ii <8> P% in the state 4>i n (pi, p 2 ) is zero. Galilean transformations, including 
global boosts and translations, are represented by unitary transformations that are local 
with respect to the interparticle tensor product structure 7ii ® H 2 and therefore do not 
affect the value of the interparticle entanglement [llj. In other words, the operator that 
performs the boost to the COM frame factors as U(-(P)) = Ui(-(P)) <g> U 2 (-(P)). 
The value (P) is invariant under the dynamics as long as the interaction is Galilean 



One boundary condition of a scattering experiment is that the in-state <pi n (formally 
the state in the limit t — > — oo) is separable with respect to the interparticle tensor 
product structure. Since (f>i n (pi,p 2 ) = 4 ) in,i{pi)4 > in,2(j>2) , one calculates that pi2(0m) = 1 
for every scattering system. The dynamics will generally not preserve this separability. 
For example, when two particle with spin scatter non-relativistically, one can show that 
even in the simplified case of central interactions and narrow momentum distributions, 
the set of S-matrices acting on the spin degrees of freedom that lead to separable out- 
states depends on the specific in-state and it is a set of lower dimension on the manifold 
of all possible symmetry-preserving and unitary S-matrices [12J. 

Another boundary condition for scattering is that the particle wave functions 
represent states that are "incoming" before the scattering. If the interaction potential 
has finite range, "incoming" suggests that the single-particle position expectation values 
in the in-state are on opposite sides of the potential region. Assuming the COM reference 
frame, the single-particle momentum expectation values are equal in magnitude and 
directed toward the potential region. Further, before the interaction the position wave 
function should have no (or essentially no) support in the potential region and the 



Pi2(4>) = / dpidp 2 dp[dp , 2 (j)(p 1 ,p2)(j)*(p[,P2)(l)(p[,P2)(l>*(PhP2) 



(3) 



invariant. 



Entanglement Mechanisms in Scattering 



5 



momentum wave functions have support only (or essentially only) on the positive semi- 
axis for one particle and on the negative semi-axis for the other. For our calculations, 
we consider the product of two Gaussian wave packets: 

(Pl-fc) 2 (P2+fc) 2 

= N 1 N 2 e ipiai e ~^T e ip2a2 e , (4) 

where iVj = (27rof )~ 1//4 , is the magnitude of the momentum of each particle in the 
COM frame, Oj are the central positions, and <7j are the momentum uncertainties for 
each particle's Gaussian. As long as a\ = —a 2 is large and when fc/cj <C 1, this wave 
function satisfies at least these heuristic notions of incoming. 

A more refined notion of incoming boundary conditions is the Hardy space 
hypothesis of A. Bohm and collaborators [13J. In that formulation, further restrictions 
are placed on the space of allowable in-states, which are defined by the preparation 
apparatus, such as an accelerator. An alternate CSCO for two particle elastic scattering 
is {P, W, S} with generalized eigenvalues of the total momentum pel, the internal 
energy w G 1R + , and relative momentum direction x = ±, respectively. Then the 
incoming wave function <f)i n (p,w,x) is a Hardy function from below in internal energy, 
i.e. it is the boundary value on the real semi-axis of a function that is analytic in the 
lower-half complex plane when w is extended to complex values. Additionally, the 
wave functions are Schwartz, giving them well-behaved smoothness and convergence 
properties in the internal energy and the total momentum. Conjugate requirements 
apply to the wave functions that represent the out-observables, which are defined by 
the detectors, but these do not enter the present analysis. It is an open question as 
to whether the requirements of the Hardy space hypothesis are consistent with the 
separability constraint on the in-state described above. However, since the Hardy- 
Schwartz spaces are dense in the Hilbert space, there will always be elements as close 
to separable as would be physically indistinguishable. 

We will not consider the intricacies of the time- dependence of the scattering 
entanglement. Instead, since the in-state particles are always unentangled, any 
entanglement in the final out-state will have been generated in the scattering event. 
The out-state (formally the state in the limit t — > +00) is found by 

4>out = S<j) in , (5) 

where 5* is the scattering operator. The exact form of the S-operator can be calculated 
for finite-range potentials by transforming to the COM-relative momentum coordinate 
system 

P = Pl+P2 (6) 

where fii = mij[m\ + m^) and solving the time independent Schrodinger equation in 
the relative momentum variable q [Hj. The S-matrix in the (p, g)-basis of the CSCO 
{P,Q} is 

(P, q\3\p', q) =S(p' - p) (t(q)6(q - q) + r(q)6(q + q')) (7) 
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The functions t(q) and r(q) are the transmission and reflection amplitudes, and unitarity 
implies |t(g)| 2 + |r(g)| 2 = 1. We also note in passing that the S-operator is a local 
operator with respect to the tensor product structure dictated by the CSCO {P, Q}, 
and so entanglement with respect to that tensor product structure is dynamically 
invariant [15] . In other words, one can define a transformed wave function <p(p, q) and 
calculate 

PpgO) = J dpdpqdp'dq'<p(p, q)(f)*(p', q)$(p' <f)<j>*(p, q'). (8) 

and one would find that p pq (4>in) = P pq {<t>out)- 

Using <£l§ and transforming back to the CSCO {Pi, A}, the out-state can be 
expressed as the sum of a transmitted and a reflected mode 



where 



and 



<t>out{Pl,P2) = 4>tra(jPuP2) + <t>ref{Pl, P2) (9) 
<f>tra{Pl,P2) = t(fl 2 Pl - //ip 2 )0m(Pl,P2) (10) 



<PrefiPl,P2) = r{fi 2 Pl ~ ^lP 2 )<fe(Pl,P 2 )- (H) 

The wave function <fe(pi,p 2 ) is the in-state wave function (/>mG°i,p 2 ) transformed by 
the reflection of the internal momentum q — > — q. One can show that 

(fa;(Pl,P2) = <j)in(Pl,P 2 ) ( 12 ) 

where (pi,p 2 ) are 

Pi = (A*i - ^2)Pi + 2//ip 2 

p 2 = 2/j, 2 pi + (/i 2 - A*i)p2- (13) 

The wave functions (ptra(Pi,P2) and re /(Pi,P 2 ) are orthogonal modes. The domain 
of support for 4>i n (pi,P2) (and therewith <ptrai.Pi ^2)) can be chosen without loss of 
generality as the region where p\ > and p 2 < 0. The domain of support of 0^G°i,p 2 ) 
(and therewith re /(pi,p 2 )) is then the region where (/ii — /i 2 )pi + 2/iip 2 > and 
2/i 2 pi + (/U 2 — Hi)p2 < 0. Remembering \i 2 = 1 — \i\ and 1 > jii > 0, one can show 
these domains have no intersection. Since the domains of support of the transmitted 
and reflected states are disjoint, the purity of the out-state is the sum of the purities of 
those two modes: 

Pl2{<Pout) = Pni&ra) + Pl 2 (0re/) (14) 

From this observation, two distinct types of entangling mechanisms can be identified. 
One source of entanglement is the superposition of the transmitted and reflected 
modes. As long as there is not perfect reflection or perfect transmission, we will find 
Pi2{4>out) < 1- The other is the entanglement within the transmitted and reflected 
modes themselves, and these mechanism can be further refined into entanglement due 
to reflection distortion and entanglement due to the variation of t(q) and r(q) with q. 
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Figure 1. These contour plots depict the bi-momentum probability densities 
\(j)f n (pi,p 2 )\ 2 (upper, left subfigure) and |0^-(p 1 ,p 2 )| 2 (the rest) for various values of 
the mass fraction /ii = m\j(vn,\ + m-i) plotted on the (pi,P2) plane. For each graph, 
the central momentum is k — 0.1 and the momentum variances are o\ = fc/10 and 
o "2 = fc/5. As the contours move outwards, each represents a reduction of probability 
density by a factor of 10. Unless the major and minor axes of the ellipses align with 
the (pi,j?2)-axes, the wave function has interparticle entanglement. 



To understand the effect of reflection, note that the transformation of momentum 
variables ffT3]) typically distorts the shape of the wave function and disrupts separability 
(see Figure 1). Because a given value of relative momentum q may overlap with a section 
in the (pi,]?2)-plane, momentum correlations ensue when the transformation q — > —q 
occurs. However, there are some special cases when this mechanism does not create 
entanglement. For the case of equal masses — — 1/2, we find (Pi,p 2 ) — > G°2,Pi) 
and so the function 4>j^(pi,P2) — * 4>in(P2,Pi) is still separable for any wave function 
4>in(P2,Pi) that satisfies the incoming boundary conditions. 

If we restrict ourselves to Gaussian in-state wave functions 4>f n {pi,P2), then we find 
that the reflected state (fi^(pi,P2) = (fifniPii^) wm i n addition be separable if 

nii/ol = m 2 /(rl, (15) 

a relationship first noted by Schulman [TJ. 

More generally, using the results for entanglement under linear transformations of 
observables found in [10], an analytic expression for the purity ^12(0^(^1,^2)) can be 
found: 

M&) = , (16) 
m - Wrt + 4^1X4/^ + (,x 2 - /iOVl) 

Using fi 2 = 1 — fii and = 0-2/0"!, this can be re-expressed as 

Pi2(#) = ° , (17) 

m y/Wi ~ I) 2 + 4/^X4(1 - tn? + (1 - 2^)2^) 




0"l/0"2 

Figure 2. This contour plot depicts the values of the function p 12 {<j&) P"7|) . The 
thick contours trace the two lines of maximum purity, corresponding to the values of 
111 and c = <J2/<J\ where the Gaussian bi-momentum wave function is separable even 
after the reflection q — » — q. Each subsequent contour represents a purity reduction of 
0.1. 

and this form highlights the fact that it is the ratio of the momenta variances, 
and not their scale, that is important. This function takes a maximum value of 
1 when either ra\ = m 2 or m\/a\ = m 2 /a 2 (see Figure 2). For these particular 
combinations of constants, the terms in the exponential that are proportional to the 
product pi x p 2 all cancel out and <j>S-(pi,P2) is again a product of Gaussians in p x and 
Pi- When mi = m 2 the variances (01,02) switch roles and <f)S-(px,p2) = (j>f n (p2,Pi), 
whereas when mi/af = m 2 /a 2 , we find the wave function is unchanged by reflection 

<©Pi,P2) = 4>f n ipi,p 2 )- 

The third mechanism for entanglement is the distortion of the wave functions of 
the two modes due to the variation of t(q) and r(q) with the relative momentum q. 
Consider the transmitted mode <f>t r a(pi,P2) — t(q)<fii n (pi,P2). Since q = p, 2 Pi ~ H1P2, 
the transmission amplitude t(q) will generally not be a separable function of p\ and 
p 2 . Each value of q corresponds to a section in the {pi-,p 2 ) plane, and these different 
sections will be given different weights when the in-state is convoluted with t(q) to 
get the wave function for the transmitted mode. Therefore, <f>tra(j>i, P2) will not be 
separable with respect to the particle momentum variables, and momentum correlations 
will ensue within the transmitted mode. A similar effect will take place in the reflected 
mode 4>ref(Pi,P 2 ) = r (<?)0i^(Pi) P2) because of the inseparability of r(q), except then 
the function <fe(pi,P2) may also be inseparable due to entanglement by the second 
mechanism described above. 
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3. Results for specific potentials 

The following subsections make explicit calculations of entanglement for specific 
potentials. In all cases, the in-state is assumed to have the form (j3J). This means that 
the in-state can be fully described by five parameters: mi, m^, k, oi, and <Ji- Instead 
of mi and 777.2, the following results will be expressed in terms of the mass fraction of 
the first particle \i\ = mi/M and the total mass variable M = mi + m 2 . 

To understand how the three different mechanisms described in the previous section 
contribute for these potentials in different parameter regions, we will consider the 
following two approximations. The first approximation is the coarse approximation, 
which we will denote as approximation (C). In (C), only the entanglement due to the 
superposition of transmission and reflection contributes. The system is effectively two 
two-level systems, and the purity takes its minimal values (corresponding to maximal 
entanglement) when the uncertainty between transmission and reflection is maximal. If 
T and R are the transmission and reflection probabilities, then 

Pu^out) =T 2 + R 2 . (18) 

Physically, the coarse approximation (C) will be a good approximation if two facts are 
true. First, the scattering amplitudes are constants over the support of the the in- 
state wave function and will be evaluated at the central momentum, i.e. t(q) — ► t(k) 
and r(q) — > r(k). The transmission probability is T = \t(k)\ 2 and the reflection 
probability is R — \r(k)\ 2 . Then the third mechanism does not apply since the scattering 
amplitudes are constants (and therefore do not disrupt separability). Second, in the 
coarse approximation (C) we neglect entanglement due to reflection distortion, which 
is reasonable if the particles have equal masses or the if the in-state is approximately 
Gaussian with masses and variances satisfying ( |T5l) . 

In the second approximation, called (C+R), we continue to assume that the 
scattering amplitudes are essentially constants, but we allow for entanglement due to 
reflection distortion. For the transmitted mode, one calculates 

Pl^tra) = \t(k)\ 4 (19) 

and for the reflected mode one calculates 

Pi2(4> ref ) = \r(k)\ 4 P i 2 (4^)- (20) 
For Gaussian in-states like (j4]), we find 

Pi2^out) = \t(k)\ 4 + \r(k)\ 4 pi 2 (^) (21) 

where P12 (</£;) is the explicit function of Oi, 0%, and /x 2 in (fTE]) and does not depend 
on the exact nature of the potential (although t(k) and r(k) do). Because Pizitffa) < 1, 
approximation (C+R) is always less than the coarsest approximation (C); reflection 
distortion can only increase entanglement in this approximation. 

For the exact results for puifiout), the four- dimensional integral (j3J) was calculated 
numerically. 
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3.1. Hard core potential 

This is the simplest case of potential scattering. In the relative variable x = X\ — £2, 
the potential has the form 



Solving the Schrodinger equation in the relative momentum gives the trivial answer 
r(q) = — 1 and t(q) = for all q. For this potential, with equal masses, the coarsest 
approximation (C) would suggest that a long time after hard core scattering there is 
no entanglement, and this result is reported in pTJ. However, as noted in [3], [8] and 
confirmed by these calculations, reflection can cause distortion. In fact, for this case the 
approximation (C+R) is exact 



Figures 1 and 2 can therefore also be considered as depicting the entangling effects of 
hard core scattering. Note that this entanglement is not dependent on any absolute 
scale such as the energy, momentum, or mass, but depends on the ratios of masses 
and variances. This kind of entanglement due to reflection (which is a factor in 
other potential scattering results below) does not disappear in the narrow wave packet 
approximation, even though the scattering amplitudes are slowing varying or not varying 
at all. 

This kind of entanglement due to momentum distribution correlations within a 
single mode can be expected to be difficult to measure compared to entanglement 
between transmission and reflection. For direct measurement, one would need a device 
to measure the momentum distribution of the scattered particles that has at least 
a resolution smaller than the in-state momentum variances. One possible scheme, 
developed for the study of atom-photon wave function entanglement in spontaneous 
emission [16], involves comparing wave function variances found in both single-particle 
and two-particle coincidence measurements of position. The Fourier transformation 
between position and momentum wave functions is local with respect to the interparticle 
tensor product structure, so measuring position entanglement is equivalent to measuring 
momentum entanglement. The practicality of this scheme would depend strongly on the 
specific nature of the system under investigation, and warrants further study. 

3.2. Dirac delta potential 
This potential has the form 




(22) 




(23) 




(24) 



The reflection amplitude is 



r{q) 



1 



1 



(25) 



kh? 



<xM>i(l-,ui) 



— I 
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Figure 3. The data points on the main plot depict the function Pi2{4>aut) an( i the 
approximation (C) and (C+R) for the case m-i = Ami for the Dirac delta potential. On 
the horizontal axis, the central momentum of the wave function k is measured in units 
of a(mi +TO2)/fi 2 . The black line represents both the coarse approximation (C) (which 
is independent of o\j '02) and the approximation (C+R) in the case a\/u2 = 1/2 (which 
satisfies the Schulman condition (fT5|) ). The dashed line is (C+R) for ci/o^ = 1 and 
the dot-dashed line is (C+R) for <J\ju2 = 2. The diamonds, circles and squares are the 
exact result computed numerically for (di = /c/10,O2 = k/5), (<Ji = k/h 1 <72 = k/5), 
and (o"i = fc/5,<72 = k/lQ), respectively. 



and the transmission amplitude is t(q) = 1 + r(g), where [i is the reduced mass 
fi = /ii/i2M. This potential has no resonances, and since the single bound state that 
occurs when a < does not participate in the elastic scattering, we can replace a — > \a\ 
without affecting any scattering entanglement results. This potential is considered in 
[21 [6], [9] for the equal mass case. 

To highlight the distinction between the exact result for the purity, approximation 
(C) and approximation (C+R), consider Figure 3 which depicts to particles scattering 
via the Dirac delta interaction when m,2 = 4m\. The figure reveals that when the 
particles have different masses, the relative momentum variance <J\j<Ji dramatically 
effects the entanglement. Unless the condition (Tl5|) is fulfilled, the entanglement is 
enhanced by the reflection mechanism for for low k because reflection dominates the 
scattering (see Figure 3 inset showing reflection probability \r(k)\ 2 ). Generally, the 
maximum entanglement (minimum purity) occurs at a value of k where the transmission 
and reflection probabilities are equal, but the location of the extreme shifts to lower k 
due to this reflection distortion, which for a given mass ratio, depends only a\jo2- 

Also note that in Figure 3, the approximation (C+R) and the exact result are very 
close, even for relatively wide Gaussian wave functions Oijk = 1/5. The variation of 
t(q) and r(q) does not lead to much wave function distortion. In the next example, 
because of resonances the scattering amplitudes vary with k at a faster rate and this 
will no longer always be true. 
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Pl2 




Figure 4. The data points on the main plot depict the function Pi2{4^ u t) an d 
approximations (C) and (C+R) for the case mi = 4mi for the double Dirac delta 
potential. On the horizontal axis, the central momentum of the wave function k is 
measured in units of a{m\ + m^jf?. The black line is approximation (C) and the 
dashed line is approximation (C+R) for crxfai — 2. The circles are the exact result 
computed numerically for a\ = fc/5. 



3.3. Double Dirac delta potential 
This potential has the form 

V(x) = a(S(x + a) + 6(x-a)). (26) 
The transmission amplitude is 

q 2 /b 2 

^ = (e 4 -<?- l) + 2iq/b + q 2 /b 2 ^ 

and r(q) = t(q) — 1 where b = (mi + m2)a/h 2 . This potential has resonant transmission 
|t(g)| 2 = 1 for particular values of the relative momentum q. For plots below, we choose 
a = 10b~ l and measure q in units of b. This case is considered in [5J for equal mass 
particles, and our methods clarify those results about the relative widths of resonances 
and the wave functions. 

In Figure 4, we contrast the two approximations and the numerically-evaluated 
exact results for m 2 = 4mi, o\ = fc/5, and o~2 = k/10. Figure 5 looks at this region 
around the first, narrowest resonance in more detail for a variety of absolute scales for 
the variances, but the same fixed ratio of variances o"i/cr 2 = 2. At resonance, the purity 
of both approximations becomes unity because only the transmission mode contributes. 
We see that in the exact results, the rapid variations of the out-state entanglement as 
found in the approximations are somewhat smoothed over. The wider the momentum 
variance in the in-state, the more pronounced this effect is. 

In agreement with [5], we also see a slight enhancement of entanglement (reduction 
of purity) for wave functions wider than the narrow resonance. As the variance scale 
get smaller, this effect disappears, and the exact results become closer and closer to 
the approximation (C+R). Note that even as the wave functions become narrower, the 
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Figure 5. This depicts an enlargement of Figure 4 around the first transmission 
resonance. The circles, squares, and diamonds are the exact result computed 
numerically for <J\/(T2 — 2 with different variance scales <j\ — fc/5, o\ — fc/10, and 
a i = fc/50, respectively. 

coarsest approximation (C) consistently underestimates the entanglement; entanglement 
due to reflection is independent of the scale of the variances and only depends on the 
masses and the ratio of the variances. Note that because of convolution with a fast- 
varying scattering amplitude, some narrow wave functions will become more entangled 
than wider wave functions, so the general statement 'resonances increase entanglement' 
should be evaluated with caution. 

4. Conclusion 

In summary, by employing symmetry methods and applying the time asymmetric 
boundary conditions, the problem of scattering entanglement in one dimension can 
be analyzed by the relative importance of three different mechanisms: two-mode 
superposition, reflection distortion, and scattering amplitude distortion. 

The overall momentum dependence of the entanglement is determined on a coarse 
scale by the two- mode effect, but if one could measure the momentum distributions 
of the two particles in the out-state, then further entanglement would be detected. 
The entanglement due to reflection is intriguing because it depends on the ratio of the 
particle masses and the ratio of the momentum variances. The effect is most pronounced 
when the more massive particle has a more certain momentum. No matter how 
sharp the initial momentum distributions are, quantum correlations ensue for Gaussian 
states unless the Schulman condition ffl5|) is satisfied or the masses are equal. As the 
distribution gets narrower, however, it would also become more difficult to measure 
this entanglement. Scattering amplitudes that vary rapidly with relative momentum on 
the scale of the variances also distort the wave function in an inseparable manner, and 
in contrast, this kind of entanglement becomes less prominent for narrow momentum 
distributions. 
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More complicated scattering problems must be considered if these methods will be 
useful for practical applications to quantum information processes with cold atoms and 
solid state devices. Generally, these will require mult i- dimensional results for identical 
particles with spin (although some one- dimensional scattering may have applications; 
see p2] for examples). In addition to possible new effects, these three mechanisms 
should still be applicable to these situations. The application of symmetry methods and 
the restriction of time-asymmetric boundary conditions remain valid, and they will be 
the starting point for future generalization of this work. 
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